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BLIND PROCESS AND RECEIVER TO DETERMINE SPACE- 
TIME PARAMETERS OF A PROPAGATION CHANNEL 

FIELD OF THE INVENTION 

This invention relates to a process and a receiver for blind 
identification of characteristic parameters of a propagation channel, 
particularly from second order statistics of the signal received by a 
sensor. 

For example, the identified parameters are the delay and 
attenuation of channel paths. 

The invention is particularly applicable in the field of mobile 
radiocommunications or ionospheric HF links. 

BACKGROUND 
Specular channel 

HF mobile or ionospheric radiocommunications are affected by 
the multiple path phenomenon. In urban areas, the transmitted signal 
is reflected and diffracted on fixed or moving obstacles present in the 
environment. For HF transmissions, the reflections are made on 
different layers of the ionosphere. In the case of propagation channels 
that can be qualified as specular, in other words the transmission 
takes place along a limited number of discrete or temporary paths 
characterized by a delay and a complex attenuation. Assuming an 
observation time compatible with the stationarity duration of the 
channel, multi-sensor reception and specular propagation, the 
expression for the multi-sensor channel is: 

c(t)=2X£(t-x k ) (1) 

where k is the index of a path, a k is the vector for which the 
components are complex attenuations of the path k for the different 



channels, t k is the delay associated with the kth path and d is the 
number of paths in a channel. 

Furthermore, if each path k is incident on the network following a 
reduced spatial diffusion cone, the expression (1) is in the following 
form 

c(t)=Z p k a(6 k )5(t-t k ) (2) 

where a(9 k ) is the input vector for an angle associated with the kth 
path and p k is the complex attenuation of the path. 

In the first case (Eq. 1), the channel is not defined by 
parameters related to the direction from which the paths arrive, each 
propagation path is defined by parameters consisting of an arrival time 
x and an « antenna vector » a. Thus, the calibration is no longer 
necessary and algorithms are no longer limited by spatial dispersion or 
coherent paths. 

In the second case (Eq. 2), paths are defined by parameters 
consisting of their directions of arrival, which assumes that the type of 
the antenna is known and therefore generally involves setting up a 
calibration to estimate the values of G k . 
Blind identification 

In active or driven systems, the channel parameters are 
calculated during a learning phase in which the transmitter transmits a 
sequence known to the receiver. 

If the propagation channel fluctuates in time, particularly due to 
movements of mobile stations or ionospheric layers, the sequence 
must be sent periodically in order to update the value of the 
parameters. 

If this type of system is efficient, the regular transmission of a 
learning sequence can cause a significant reduction in the effective 
throughput. 



For example, in the STANAG 4285 standard for cooperative 
Ionospheric HF transmissions, half of the transmitted symbols are 
learning symbols. 

Prior art also includes different blind methods and systems, in 
which parameters are estimated starting from statistics of the received 
signal, without any advance knowledge of the learning sequence. 
Second order techniques 

For example, the proposed techniques simply use second order 
statistics (space-time covariance matrix) of the received signal. 
Second order algorithms have better convergence properties than 
higher order algorithms, i.e. the variance of second order estimators 
for a given number of symbols, is usually less than the variance of 
estimators with higher orders. Furthermore, they have fewer local 
optimisation problems than techniques with higher orders. 

Various methods and algorithms have been developed, such as 
those described in document entitled "Multichannel Blind Identification 
: from subspace to maximum likelihood methods" by L.Tong and 
S.Perreau; Proceedings of the IEEE, 86(10): 1951-1967, October 
1998. One of the disadvantages of these algorithms is that the length 
of the global transmission channel has to be known, and most of them 
are intolerant to an error in the estimate of this length. When the 
bandwidth of the transmission is limited, the channel length is only 
defined approximately and these algorithms can no longer be used. 
State of the art for blind parametric identification 

Most work on blind identification of characteristic propagation 
parameters uses decoupled algorithms, i.e. algorithms that 
independently estimate arrival times and antenna vectors or arrival 
directions. These include direction finding algorithms to estimate 
directions of arrival, and coherent source separation algorithms to 
estimate antenna vectors. A lot of work has also be done on 
estimating arrival times. 



Joint estimating methods can improve the precision and 
resolution of estimators, such that parameters can be estimated even 
when delays (or angles) are very similar. 

The main work done on blind joint estimating of parameters 
(0,x) is described in the following references: 

• "Identification spatio-temporelle de canaux de propagation a trajets 
multiples" ("Space-time identification of multipath propagation 
channels) by J.Gouffraud, PhD thesis, Ecole Normale Superieure de 
Cachan, 1997, and 

• "Improved blind channel identification using a parametric approach" 
by M.C.Vanderveen and A.Paulraj, IEEE Communications Letters, 
pages 226-228, August 1998. 

The idea is that criteria used to make a blind estimate of the 
pulse response may be directly minimized as a function of the angles 
and delays, using subspace type criteria like those described in the 
"Subspace methods for the blind identification of multichannel FIR 
filters" paper by E. Moulines, P. Duhamel, J-F. Cardoso and S. 
Mayrargue, IEEE Trans, on signal Processing, 43(2): 516-525, 
February 1995. These algorithms require that the transmis- 
sion/reception filter is known in advance and that the antenna is 
calibrated. 

Furthermore, an example of a method of making a joint blind 
estimate of parameters (a,x) is described in the "Methods for blind 
equalization and resolution of overlapping echoes of unknown shape" 
paper by A.Swindlehurst and J.Gunther, IEEE Trans, on Signal 
Processing, 47(5): 1245-1254, May 1999. The authors work in the 
frequency range and propose an iterative IQML type algorithm to give 
an approximate solution with a maximum probability and an explicit 
initialisation algorithm based on the ESPRIT algorithm. These algo- 
rithms do not require advance knowledge of the transmission filter, but 



do have the disadvantage that a Fourier transformation is necessary 
before processing. 

SUMMARY OF THE INVENTION 

The invention relates to a blind or partially blind process that 
can determine at least the characteristic space-time parameters of a 
propagation channel in a system comprising at least one sensor for 
reception of signals y(t), characterized in that it comprises at least one 
step in which the specular type structure of the channel is taken into 
account and a step for the joint determination of space-time 
parameters such as antenna vectors (a) and/or time vectors (x) 
starting from second order statistics of the signals received on one or 
several sensors. 

According to one embodiment, it comprises a step in which the 
received signal y(t) is oversampled, for example with a sampling 
period equal to T/p, where T is the symbol period. 

According to another embodiment, the process comprises a 
reception step on at least two sensors and an oversampling step. 

For example, the process according to the invention may be 
applied for standard links for equalization, positioning and spatial 
filtering purposes, or for monitoring the spectrum of a propagation 
channel for positioning purposes starting from one or several HF 
stations. 

In particular, the invention has the following advantages: 

• there is no need to transmit a learning sequence that may reduce 
the effective transmission throughput, or to know the length of the 
global transmission channel, 

• it is tolerant to an overestimate of the order of the transmission 
channel. 



BRIEF DESCRIPTION OF THE DRAWINGS 

Other advantages and characteristics of the invention will 
become clear after reading the following description with reference to 
the single figure given for illustrative purposes and in no way limitative, 
and that shows the schematic for a radio-communication receiver 
according to the invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 

We will mention a few remainders, assumptions and definitions 

before going on to describe the steps in the process according to the 

invention. 

Assumptions and model 

The single sensor signal model 

The transmitted signal is in the form 

x(t)=Z g(t-IT)s(l) (3) 

/eZ 

where g(t) is the transmission/reception filter and values in {s(l)} are 
information symbols that were sent at the symbol period T. 

The system assumes that the carrier frequency and the symbol 
period are estimated in advance. The signal received in base band is 
written 

Y(t)=Z h(t-IT)s(l)+w(t) (4) 

/eZ 

where the values of {s(l)} are information symbols that were sent at the 
symbol period T, h(t) is the pulse response of the transmission 
channel, and w(t) is an additive noise. The received signal sampled at 
the symbol period follows the discrete model: 

y(k) = y(kT)Z h(k-l)s(l)+w(k) (5) 

leZ 



where h(l) = h(IT) (assuming that the sampling phase is equal to 0). 
The multichannel system 

Figure 1 represents a multichannel reception system 
comprising a number q of reception sensors and a means of 
oversampling the signal s(k) corresponding to path index k. 

Obviously, without going outside the framework of the invention, 
the process can also apply to receivers comprising several sensors 
without any means of oversampling, or comprising a single sensor and 
a means of oversampling the received signal. 

In the example embodiment given below, the device considered 
is a multichannel system combining spatial diversity and time diversity, 
i.e the signal is received on a network of q sensors Ci each having a 
response h (i )(t) and it is then oversampled with respect to the symbol 
period T, at the rate T/p as shown in figure 1 . Oversampling is done 
using an appropriate device. The global pulse response of the 
transmission channel and the multi-sensor observations vector are 
defined by the following qx1 column vectors: 
h(t) = [h (1) (t),..„ h (q) (t)] T (6) 

y(t) = [y<i)(t),..,y(q)(t)] T (7) 

The subscript indicates the sensor number and the superscript 
T represents the transpose of a matrix. Bold characters indicate a 
spatial dimension. The samples obtained by sampling at rate T/p are 
vectorized by collecting signals sampled in the same sampling phase 
within a symbol period. More precisely, a path is associated with each 
of the p different sampling phases at the symbol rate, thus creating a 
multipath channel with p outputs on each sensor. This gives a 
multipath channel with p.q outputs. The global pulse response of the 
transmission channel and the multipath observations vector are 
defined by the following column vectors p.q x1: 

h(k) = [h (1) (kT)... h (p) (kT)] T (8) 

y(k) = [y (1) (kT)... y (p) (kT)] T (9) 



The superscript indicates the sampling phase, 

h {i) (kT) = h(kT +(i - 1 )T/p) f y (i) (kT) = y(kT + (i - 1 )T/p). 
Assumptions 

The process according to the invention is based on a number of 
assumptions, and particularly the following. 
Symbols source 

H 1 - For example, process { s(k) } is a complex white noise that may 
also be circular at the second order of unit variance, 

E [s(k)] = 0, E [s(k)sG) H ] = 5(k - j), E [s(k)sG)] = 0 (10) 
Noise 

H 2 - Noise w(k) is a Gaussian, circular, white process stationary in 
space and time with variance a 2 , 

E [w(k)w(k) H ] = a 2 Ipq (11) 
where Ipq is the identity matrix 
Transmission/reception filter 

H 3 - The transmission/reception filter g(t) has a finite support, g(t) = 0 
for t e [0, LgT[ where Lg corresponds to the filter length and T to the 
sampling period. 

H 4- The transmission/reception filter g(t) <-> G(f) is approximately as 
follows for a limited band, with band B g 

6(f) *0 Vf £ [-B gi B g ] ( 12 ) 
and G(f) does not cancel out within a given band (3, 

G(f)*0 Vf e [-B g ,B g ] (13) 
Propagation channel 

H 5 - The signal is a narrow band for the sensor network, i.e. it is 
assumed that the propagation time for the signal to pass from one 
sensor to the next is very much less than the inverse of the signal 
band. 

H6 - The channel is constant over an observation interval. The 
Doppler shift and carrier residues are assumed to be negligible. 



H 7- The time difference between the first and last path is finite, and it 
has an upper limit: Ar max . 

where d is the number of paths and r k is the delay and a k is the 

vector of q x 1 dimension antennas associated with the kth path, the 
channel response being expressed as follows. 

h(t)=£ a k g(t-T k ) (14) 



By sampling at T/p, 

h r fz> ±h(kT/p)z- k =ta k g T (T k ;z) (15) 

- JU-oo k=I - 

P P 

where 

Q T (T;z)AZg(n--T)z- n (16) 
7 P 

is the transform into z of the transmission/reception filter shifted by 

r and sampled at T/p. 
Furthermore, h(z) is the response in z to the symbol rate of the 
combined system, 

h(z)= ± h(kT)z k (17) 

k=-<£> 

The received signal follows the equivalent single source multi-output 
discrete model: 

y(k) = [h(z)]s(k) + w(k) (18) 
it is interesting to note that 

h I (z)=2 h»(zV (19) 

where h {,) (z) = ^r=-« h (1) (k)z k are the multi-phase components of the 
oversampled channel. 



According to assumptions H3 and H7, the pulse response of the 
global multi-sensor transmission channel has a finite length. 
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Therefore, the degree of the discrete responses h(z) and h T / P (z) is 
finite. We note L = deg(h) i.e. the degree of at least one of the paths 
is equal to L, the degrees of the others may be lower, 

h(z) = h(0) + h(1)z* 1 + ... + h(L)z L 
Similarly, we denote L T / P = deg(h T / P ). We have pL < L T / P < p(L + 1)-1. 
In order to simplify the presentation, we assume that 

L T / P +1=p(L+1), 

Furthermore, second order algorithms are based on the following 
diversity condition for the global transmission channel: 
H 8 The polynomial vector h(z) is irreducible, i.e. these components do 
not have any common roots. 
Some notations and tools 

If {b (1) (z),..., b (r) (z)} are r scalar polynomials, 

b(z) = [b (1) (z), b (r) (z)] T is a polynomial vector (a vector for which 
the coordinates are polynomials) with dimension r. The maximum 
degree of the polynomials is called the order or the degree of the 
vector deg(b(z)) 

{b (1) (z), b (r) (z)}. If L b is the degree of b(z), then 



where b(k) = [b (1) (k) b (r) (k)] T . The vector associated with b(z), i.e. 
containing all its coefficients, is denoted by an arrow, 



with dimension r(L b + 1) x 1. Therefore, in particular h 
corresponds to the dimension vector 
q.p(L + 1) x1 associated with h(z) and h T / P (z) 
i.e. containing all their coefficients, 



b(z) = X b(k)z k 



(20) 



£=[b(0) T ...b(L b ) T ] T 



(21) 



h =[h(0) T ...h(L) T ] T 



(22) 
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= [h r (0) T ...h r (L r ) T ] T 



(23) 



Furthermore, let K be an integer defining the number of 
multidimensional observations, K (b) is the Sylvester matrix (i.e. 
Toeplitz per block) with dimension r(K + 1 ) x (K + U + 1 ) 

associated with the polynomial b(z), and defined by 



T k (b) = 



'b(O) b(L b ) 0 
0 b(0) b(L b ) 



0 



0 b(0) 



0 



(24) 



Finally, let B(z) be a polynomial matrix with dimension r x D and 
degree L b , 

B(z) = [bi(z),..., b D (z)] 
where bj(z) is a polynomial vector with dimension r. We will use 

B = [5i b D ] 

to denote the associated matrix with dimensions r(U + 1) x D, and 
Dk(B) to denote the Sylvester block matrix with dimensions r(K + 1 ) x 
(K + L b + 1)D defined by 



D k (B) = [T k (bi) ) ... I T k (b D )] 



(25) 



^min(A) denotes the smallest eigenvalue of matrix A and vp m i n (A) is the 
eigenvector associated with X m \„(A). 
Process according to the invention 

The process according to the invention is based particularly on 
taking account of the specular type channel structure to determine the 
specific parameters of a signal propagation or transmission channel, 
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particularly in a communication system like that shown in figure 1. For 
example, the parameters determined jointly may be the time vector (x) 
of delay vectors (tk), and the antenna vector (a) of vectors (a k ). 

A first alternative embodiment of the process or the "a priori 
process" is particularly suitable when the transmission/reception filter 
is known. 

A second alternative embodiment of the process or "no a priori 
process" is more suitable for the case in which the characteristics of 
the transmission/reception filter are not known. 
A priori process 

The process assumes that the transmission filter g(t) is known. 
The various steps in the a priori process or algorithm according to the 
invention are based on the fact that equation (15) is written in 
"vectorial" form as follows: 



h is the vector with dimension q.p.(L + 1) containing all coefficients of 
the response hj/p(z) } a is the vector with dimension q.d containing all 



and Gu.dM is a matrix with dimension q.p(L + 1) x q.d containing 
delayed versions sampled at T/p of the transmission filter: 



h = G L|d (T)a 



(26) 



antenna vectors a = [a^—a]] 7 



T 



g(0-Tj 

T 



GL,d(t) = 



p 



p 



®lq 



(27) 



g(L + l)T---r } 



g(L + I)T--~ 



P 



P 
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l q is the identity matrix with dimension q. We will use x to denote the 
vector with dimension d containing all delays x = [r,,...,rj 

For example, the different steps in the process could be as follows: 

1 . Choose the length of the global propagation channel taking account 
of transmission and reception, L>L g +Ar max and the value of the 
number of observations K > L where Ar max is the largest possible 
value of the relative delay between two paths. 

2. Calculate the sub-space criterion matrix Q L (R), for example as 
follows: 

- Calculate the empirical estimate of the covariance matrix: 

R=^Zy«m«(k) H (28) 
N-K So 

where 

y K (k) = [y(0) T y(T/p) T . . . y(KT - T/p) T ] T (29) 

= [y(0) T ...y(K) T ]J g(K + 1) (30) 

and let N-K be the number of independent observations. 

- Build the matrix G L for example using the eigenvectors associated 
with the pq(K + 1) - (K + L + 1) smallest eigenvalues of R. Gj- can 
also be obtained from the breakdown of the data matrix 

Y±\y(K)---y(N )].) into eigenelements. Then form the projection 
matrix fl £ = G L G- L 

- Form the sub-space criterion matrix Q Z (R) - D i <Jl L )D- L (f{ i ) H 
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3. Estimate the number of paths, for example d , by using the method 
described in the "Detection of signals by information theory criteria" 
reference by M. Wax and T. Kailath, IEEE Trans.on Acoust. and 
Sig. Proa, 33(2) : 387-392, April 1985. 

4. Estimate the delays 



f =argmin J^^(r) 



where 



. f , KJG U (r) H Q L (R)G u (r)) 
UiT) ~ KJG u (r)"G u (r)) 



This is done by applying the iterative process described below: 

- Let d = 1. Obtain an estimate f, of r 01 by minimizing J Li {r x ) 

- Let d = 2. Obtain and estimate f 2 of r 02 by minimizing J u ({r,,r 2 ]) 

as a function of r 2 and keeping r, equal to the estimate f, obtained in the 
first step. Then, obtain a new estimate of [r 0 ,,r 02 ] using a 
multidimensional minimization of J Lji (\t l ,r 2 \) 
and using the values [f p f 2 ] that were obtained above as initial values. 

- Iterate this process until d = d , for all possible paths in order to obtain 
the delay vector. 

5. Estimate the value of the antenna vector by 

a = vp min (G u (r) w Q L (R)G u (f)) 

6. Create the estimate of the pulse response: h = G Lk (i)a 
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No a priori process 

The second alternative of the process corresponds to a 
parameter estimating algorithm that is entirely or mostly blind. As in 
the first alternative embodiment, the associated algorithm takes 
account of the specular nature of the propagation channel in order to 
determine its parameters, and particularly the antenna vector and the 
characteristics time vector. Unlike the first variant, it does not require 
advanced knowledge of the transmission filter. It simply assumes that 
the transmission filter has a limited band and that its band is known 
approximately. 

For example, the various steps in the process or algorithm are as 
follows: 

1. Estimate Lg , the length of the transmission filter, and choose L> 
L g + Aimax and K > L 

2. Apply the sub-space method described in item 2 of the a priori 
process: 

- Estimate the covariance matrix R , for example using the steps 
described above corresponding to equation (Eq. 28). 

- Calculate the projection matrix on the noise space f[ £ starting from 

the eigenvectors associated with the pq(K + 1)-(K+Z+1) null 
eigenvalues of the matrix R . 

- Form the matrix Q t (R)= D L (f[ L )D L (fl L ) H 

- Obtain h the eigenvector associated with the smallest eigenvalue of 
the matrix Q Z (R) 

and form h T (z) 
p 

3. Form the new parametric criterion: 

- Choose the value of the channel length R >L T = p(Z + 1) - 1 
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and form the matrix tr (h T ) 



p 



- Calculate the matrix Q ; containing eigenvectors associated with the 

T 



q(R + 1 ) - (R + L T + 1 ) smallest eigenvalues of t r ( h T ) 



H 



P 



P 



- Form the matrix D , ( Q , ) 

L r L T 



4. Estimate the number of paths d , for example using the method 
referenced above for step 3 in the first alternative embodiment. 

5. Determine the minimum value of the band of the transmission filter: 
P, usually (3 = 1/T. 

6. Choose v (t), a continuous filter with band limited to the p band. 
Form the filter v(t), 




0<t<Lv<Lg 
elsewhere 



(32) 



7. Estimate the delays 



r = argminl - .( r) 



where 



KJV L 3 (x) H D iL (Q L T )D L r_(Q L T )"V u (t)) 
h j (V : L ' (33) 
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using the iterative process proposed for the parametric subspace 
method. (The matrix V z ^(r) is defined in the same way as the matrix 

G u (t) but from the filter v(t).) 

8. Estimate antenna vectors: 

a = vP^ (V u (r) H D iL (Q z r )D lL (Q l t) h V u (f)) (34) 
p p p p 

Possible applications 

The process according to the invention is applied, for example 
in the domains mentioned below for illustrative purposes and in no 
way limitative. 
Cellular networks 

For cellular transmissions, the transmission standard is known. 
Therefore, the standard transmission filter is known and it is preferable 
to apply the a priori process using advance knowledge of the 
transmission/reception filter. The potential applications are as follows: 
Equalization 

Propagation acts as time convolution. Except for very low 
throughput communications (a few hundred bits per second), this time 
dispersion induced by propagation creates interference between 
symbols (IES). In order to best recover the transmitted symbols, the 
dispersive effect of the channel must be compensated by an 
equalization step like that described in the book entitled "Digital 
Communications", written by J.G. Proakis and published by McGraw 
Hill editions, 1989. If the transmission/reception filter is known, the 
global response can be built up starting from knowledge of 
propagation channel parameters, and an equalization step can be 
applied. The advantage of this approach is that in general, by 
estimating propagation parameters rather than the global response, 
the number of parameters to be estimated is smaller so that it is 
usually possible to reduce the variance of the estimators. The 
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identification quality and equalization performances can thus be 
Improved significantly. Furthermore, some characteristic parameters 
of the propagation channel (arrival delays and angles) do not change 
very much from one frame to the next. Parametric modelling makes it 
more efficient to monitor variations of channel parameters and further 
improves the identification quality, for example such as described in 
the "Multi-channel MLSE equalizer with parametric FIR channel 
identification" reference by J.T. Chen and A. Paulraj, in VTC Proa, 
pages 710-713, 1997. 
Positioning 

Positioning of the transmitter is a classical problem in military 
applications, and has recently become the subject of close attention 
for civil applications in mobile communications. Positioning of mobiles 
in a cellular network makes it possible to provide many services such 
as positioning of emergency calls, management of charging, detection 
of fraud, assistance with cellular planning, improvement in organizing 
the "handover" as described in " Position location using wireless 
communications on highways of the future" reference by 
T.S.Rappaport, J.H.Reed and B.D.Woerner, IEEE Commun. Mag., 
October 1996. 

Several techniques can be used for positioning. The most 
classical is triangularization. In this case, positioning is done using 
two or three separate receivers (base stations) each estimating the 
angles and/or arrival times of the transmitted signal as described in the 
"Mobile Localization in a GSM Network "reference by H.EI Nahas, PhD 
thesis ENST, June 1999. This system has some disadvantages, 
which is why positioning algorithms from a single station were 
developed. These positioning techniques always require prior 
knowledge of propagation characteristics. For example, a positioning 
process is described in the reference entitled "Procedes de 
localisation spatio-temporel de radiomobiles en milieu urbain " (Space- 
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time positioning processes for radiomobiles in urban environment) by 
M. Chenu-Tournier, A. Ferreol and J-J. Monot, Technical report, TCC, 
1997. 

Antenna filtering - SDMA 

Knowledge of spatial characteristics of the different paths may 
be used to focalise transmission in the main arrival direction and thus 
minimize the effect of multiple paths on receivers (uplink bearnformer 
design, downlink transmission design). 
Spectrum monitoring 

In the context of a spectrum radio-electrical analysis, the 
receiver has no prior knowledge about the transmitter. The reception 
system scans the spectrum, intercepts a signal transmitted by an 
unknown transmitter and analyses it to extract some information from 
it to be able to listen to the transmission or determine the position of 
the transmitter. In this context, it is unlikely that the transmission filter 
will be known and the no a priori process is the most suitable for 
estimating propagation channel parameters (x ; a). In this case, since 
the transmission filter is not known, it is impossible to equalize the 
received signal, but knowledge of parameters (x ; a) would make it 
possible to use the applications mentioned below. 
Positioning 

The main application is positioning. In HF, positioning with a 
single station conventionally uses an initial step of ionospheric 
forecasts. Knowing the angle of incidence and the electronic density 
profile (by ionospheric forecasting), a ray plotting method using 
Descartes law can be used to reconstitute the path of the wave 
received on the sensors network according to one of the following 
references: 

• "Comparison of the fixing accuracy of single-station location of long 
range transmitters, by H.C.Horing , IEEE Proceedings, 37(3) : 173- 
176, June 1990. 
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• "loniospheric modelling in support of single station location of long 
range transmitters by LF.McNamara, Journal of Atmospheric and 
terrestrial Physics, 50(9) : 781-795, 1988. 

Knowledge of arrival times can also eliminate some ambiguities due to 

multi-skip. 

Others 

Spatial filtering is possible if antenna vectors are known. 
Knowledge of delays is a means of recombining propagation paths in 
phase. The use of one or both of the two known values may be a first 
step in implementing a signal analysis system. 



